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ABSTRACT 



Context. The role of magnetic fields during the protostellar phase of high-mass star-formation is a debated topic. In particular, it is still 
unclear how magnetic fields influence the formation and dynamic of disks and outflows. Most current information on magnetic fields 
close to high-mass protostars comes from H2O and OH maser observations. Recently, the first 6.7 GHz methanol maser polarization 
observations were made, and they reveal strong and ordered magnetic fields. 

Aims. The morphology of the magnetic field during high-mass star-formation needs to be investigated on small scales, which can only 
be done using very long baseline interferometry observations. The massive star-forming region W75N contains three radio sources and 
associated masers, while a large-scale molecular bipolar outflow is also present. Polarization observations of the 6.7 GHz methanol 
masers at high angular resolution probe the strength and structure of the magnetic field and determine its relation to the outflow. 
Methods. Eight of the European VLBI network antennas were used to measure the linear polarization and Zeeman-splitting of the 
6.7 GHz methanol masers in the star-forming region W75N. 

Results. We detected 10 methanol maser features, 4 of which were undetected in previous work. All arise near the source VLA 1 
of W75N. The linear polarization of the masers reveals a tightly ordered magnetic field over more than 2000 AU around VLA 1 
that is exactly aligned with the large-scale molecular outflow. This is consistent with the twisted magnetic field model proposed 
for explaining dust polarization observations. The Zeeman-splitting measured on 3 of the maser features indicates a dynamically 
important magnetic field in the maser region of the order of 50 mG. We suggest VLA 1 is the powering sources of the bipolar outflow. 
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1. Introduction 

Magnetic fields are attributed an important role in star forma- 
tion, in particular in halting the collapse, transferring angular 
momentum, and powering the outflows. Gravitational collapse 
of molecular clouds proceeds primarily along the magnetic 
field lines, giving rise to large rotating disk or torus structures 
orthogonal to the magnetic field (Matsumoto & Tomisaka 
1200 4). Consequently, the molecular bipolar outflows, which 
originate in the protostar, are driven parallel to the magnetic 
field (e.g., Cohen 2005, McKee & Ostriker |2007] l. Despite it is 
largely accepted for low-mass star-formation, this picture is still 
under debate for the formation of high-mass stars. 

The best probes of the magnetic field on the smallest 
scales in the high-mass star-forming environment are maser 
lines. Their bright and narrow spectral line emission is ideal 
for detecting both the Zeeman-splitting and the direction of 
the magnetic field. So far, the most investigated sources are 
H2O and OH masers, revealing an ordered structure in the 
magnetic field and field strengths between 10 and 600 mG 
(e.g., Vlemmings et al. I2006 a^ and a few mG (e.g.. Fish et al. 
I2OO5I 1. respectively. Even though methanol (CH3OH) is the 
most abundant of the massive star-formation maser species (e.g. 
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Pestalozzi l2007b . the first 6.7 GHz (5i 60A+), methanol 
maser linear polarization observations have only recently been 
made (e.g., Vlemmings |2006bt Green et al. '2007; Dodson 
2008). Like H2O, CH3OH is a non-paramagnetic molecule, 
and thus any linear polarization is a few percent (~2-3%), 
while circular polarization is well below 1%. Still, Vlemmings 
(2008) detected 6.7 GHz methanol maser Zeeman-splitting in 
a sample of 24 massive star-forming regions, including W75N, 
for the first time using the Effelsberg telescope. This revealed 
a magnetic field strength in the methanol region of about 20 mG. 

In Fig. 1 we show the area of interest of this letter. VLA 1 and 
VLA 2, which were studied in detail by Torrelles et al. ( 119971 ). 
are part of the region W75N(B) in the active high-mass star- 
forming region W75N (Hunter et al. 119941 ) located at a dis- 
tance of 2 kpc. A large-scale, high-velocity outflow, with an 
extension greater than 3 pc and a total molecular mass greater 
than 255 Mq, was detected and mapped at different frequencies 
(e.g.. Shepherd et al. 120031 and references therein). Shepherd et 
al. (2003 ) propose a multi-outflows scenario where VLA 2 may 
drive the large-scale outflow, and VLA 1 is the centre of another 
small flow. However, so far, it has been impossible to determine 
which is the main powering source of the 3 pc outflow. 

Several maser species were detected in W75N(B): H2O 
(e.g., Torre lles et al. [T997] ). OH (e.g., Baart et al 1986 Slysh 
et aL l2002] ). and CH3OH (Minier et al. lIOOOl hereafter MBC). 
The methanol masers were detected in a linear structure and off- 
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Table 1. Results. 
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Features RA Dec Peak Flux S(v) V,^, P, x AVz By 

Density (I) 

(J2000) (J2000) (Jybeam-') (Jy) (kms-') (%) (°) (m.s-') (mG) 



Al 


20h38m36.412s 


42d37m35.285s 


15.54 


20 


3.54 


2.0 ±0.4 


-22 ± 1 






A2 


20h38m36.409s 


42d37m35.283s 


48.25 


72 


4.07 


4.0 ± 1.2 


-20 ±8 






A3 


20h38m36.405s 


42d37m35.236s 


67.24 


84 


4.60 


4.5 ±0.3 


-9.6 ±0.3 






A4 


20h38m36.400s 


42d37m35.165s 


47.58 


62 


5.82 


0.9 ±0.3 


-6 ±6 


0.80 ±0.03 


16.3 ±0.6 


A5 


201i38m36.401s 


42d37m35.154s 


39.39 


51 


5.12 


3.5 ± 1.5 


-2± 12 


0.75 ±0.13 


15.2 ±2.6 


Bl 


20h38m36.417s 


42d37m34.846s 


95.38 


107 


7.23 


1.5 ±0.5 


-18±9 


0.53 ± 0.04 


10.9 ±0.8 


CI 


201i38m36.418s 


42d37m34.693s 


2.33 


3 


9.34 










C2 


201i38m36.420s 


42d37m34.545s 


7.02 


8 


6.88 










C3 


201i38m36.411s 


42d37m34.376s 


4.11 


5 


9.16 


1.3 ±0.3 


-25 ±4 






C4 


20h38m36.403s 


42d37m34.323s 


10.83 


11 


9.51 


1.4 ±0.3 


-30 ±21 







set from the compact clump at northwest (group A) and south 
(group B) of VLA 1, respectively (MBC). VLA2 is the place 
where the most intensive OH flare took place (Alakov et al. 
I2005I I. Other OH maser emission sites are situated on a ring 
structure around the radio sources (Hutawarakorn et al. l2002l l. 

To investigate the magnetic field of this region, a number 
of OH maser polarization observations have been made (e.g., 
Hutawarakorn et al. l2002l Slysh et al. 2002; Fish & Reid|2nQ7ll. 
The magnetic field strength obtained from the OH maser po- 
larization was 7 mG (e.g., Hutawarakorn et al. 120021 Slysh et al. 
I2002I I. During the OH maser flare near VLA 2, a strong magnetic 
field of more than 40 mG was detected in several maser spots 
(Slysh & Migenes 2007 ). Effelsberg observations of the 6.7 GHz 
methanol masers indicate a magnetic field B\\ = 5.7 mG and 

= 9.5 mG in the A and B groups, respectively (Vlemmings 
I2U08 T). 

After investigating the linear polarization of OH masers (1.6 
and 1.7 GHz), Hutawarakorn & Cohen (1996) suggested that 
the magnetic field is oriented along the outflow. However, OH 
masers at 1 .6 and 1 .7 GHz are very susceptible to both internal 
and external Faraday rotation, making a direct determination of 
the intrinsic magnetic field structure difficult. Hence, the study 
of the linear polarization of the 6.7 GHz methanol maser can 
provide more precise polarization angles and, consequently, a 
more accurate determination of the magnetic field than obtained 
at lower frequencies. Here we present the first European VLBI 
network (EVN) 6.7 GHz methanol maser polarization observa- 
tions. 

2. Observations and data reduction 

W75N(B) was observed in full polarization spectral mode at 
6668.518 MHz with 8 of the EVN antennas on June 14* 2008. 
The bandwidth was 2 MHz, providing a velocity coverage 
of ~ 100 km s"'. The data were correlated twice. First, with 
modest spectral resolution (512 channels; channel spacing 
~ 0.2 km s"'), which enabled us to generate all 4 polarization 
combinations (RR, LL, RL, LR). The second correlation run 
with high spectral resolution (1024 channels; channel spacing 
~ 0.1 km s ') only provided the circular polarization combi- 
nations (LL, RR). This high spectral resolution dataset was 
used to attempt the detection of the Zeeman-splitting with a 
similar RCP-LCP cross-correlation method that was used for the 
Effelsberg observations by Vlemmings (|2008j). The observation 
time was 8 hours, including overheads on the calibrators 3C48, 
J2007+4029, 3C454.3, and J2202+4216. The last two were also 
used for polarization calibration. 



The data were reduced using the Astronomical Image 
Processing Software package (AIPS). The initial polarization 
calibration was performed on J2202-H4216. Fringe-fitting 
and self-calibration were performed on the brightest spectral 
channel, corresponding to the strongest maser feature (~90 
Jybeam"' at Visr = 7.23 km s"^). The observations of the polar- 
ization calibrators and the unpolarized calibrator J2007 -1-4029 
indicate the remaining polarization leakage after calibration is 
much less than 1 %. Then we imaged the I, Q, and U cubes (2048 
X 2048 mas) of the source W75N(B) using the task IMAGR 
(beam of 8.75x6.45 mas). The data cubes were combined to 
produce cubes of polarized intensity and polarization angle. 
With the exception of the bandpass table, which had to be 
remade, the calibration solutions were copied from the modest 
resolution dataset and applied to the high spectral resolution 
data. From this dataset we produced RR and LL image cubes. 

Our observations were obtained between two VLA po- 
larization calibration observation^ during which the linear 
polarization angle of J2202-t-4216 changed abruptly from -30° 
to -80°. However, we were able to calibrate the linear polar- 
ization angle using 3C454.3, which had a constant polarization 
angle of 6?9 from VLA monitoring during our observation'. 

The absolute position of the brightest feature was obtained 
through fringe rate mapping using the AIPS task FRMAP. 
We determined its position to be aaooo = 20''38™36'.417, 
^2000 - -h42°37'34''846 with an absolute position error of about 
20 mas. The relative position errors depend on the maser flux of 
each feature, but are better than 1 mas. 



3. Results 

We have detected 10 methanol maser features (see Table[T]), with 
a minimum signal-to-noise ratio of approximately 80. The peak 
flux and the flux density (S (v)) were obtained from a Gaussian 
fit to the brightest emission channel. The velocity with respect to 
the local standard of rest (Visi) was obtained from a fit to the fea- 
ture spectrum, as was the linear polarization fraction (Pi). The 
error on the linear polarization angle (x) includes the formal er- 
ror due to the thermal noise (Wardle & Kjonberg J_974). Figure 
[T] shows the positions of methanol maser features. They can be 
divided into three groups: A, B, and C, following the naming 
convention adopted by MBC. Using three EVN antennas in May 
1997, they detected 12 methanol maser features in W75N(B) di- 
vided into two groups: A, which is located at the northwest of 

' http://www.aoc.nrao.edu/astro/calib/polar 
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Fig. 1. Positions of methanol and water masers superimposed on the 1.3 cm continuum contour map of the VLAl thermal jet and 
VLA2 (Torrelles et al. 119971 ). Contours are -3, -2, 2, 3, 4, 5, 6, 7, 8 x 0.16 mJybeam"'. Blue triangles indicate the positions of the 
water masers as reported by Torrelles et al. (|1997| l. Red circles indicate the position of the 10 methanol maser features with their 
linear polarization vectors (40 mas correspond to a linear polarization fraction of 1%). The two arrows indicate the direction of the 
bipolar outflow (66°) and the parallel dashed lines the magnetic field lines (73°+ 10°) as derived from the linear polarization. 



VLA 1, and B to the south of VLA 1, with velocities between 3 
and 7 kms 

As found by MBC, the features A1-A5, with velocities 
between 3 and 6 kms"', are located along a linear structure. 
Although the absolute positions in MBC are uncertain by up to 
1 arcsec, we were able to identify 5 features detected by MBC 
as A1-A5. Their velocities are slightly higher (+0.04kms"'< 
A Visr < +0.16kms ') than those reported previously. Only fea- 
ture A5 (Visi— 5.12kms"') has a lower velocity ( 0.5 kms"'). 

Feature Bl, which is located southeast of the group A at a dis- 
tance of about 400 mas (~800AU at a distance of 2kpc), has 
a velocity of 7.23 km s"' . It can be associated with the brightest 
feature (peak flux = 71.1 Jy) of MBC. Comparing the separations 
between A1-A5 and Bl in the two epochs, we find a small dif- 
ference in their relative positions, implying a separation velocity 
of 0.93 mas yr"' ~ 9.5 kms"' . The group C, located south of 
the feature Bl, is composed of four methanol maser features un- 
detected by MBC. The features of this group (C1-C4) seem be 
along an arc structure with Visr between 6 and 10km s"'. 

Linear polarization was detected in 8 of the 10 maser fea- 
tures given in Table [T] The highest fractional linear polarization 
(Pi ~ 2 - 5%) was detected in group A. The linear polariza- 
tion fraction of the masers in group B and C is slightly lower 
(P\ ~ 1 - 2%). No dependence of Pi on maser brightness is 
found. The weighted polarization angle of methanol maser vec- 



tors is (x) - -17° ± 10°. Assuming the interstellar electron den- 
sity and the magnetic field are ~ 0.012 cm"^ and By =s 2yuG, 
respectively, the Faraday rotation at 6.7 GHz is w 5°. Zeeman- 
splitting was detected in 3 of 10 maser features (A4, A5, andBl). 
Figure |2] shows the total intensity and the circular polarization 
spectrum for the brightest feature. Columns 9 and 10 of Table 
[1] contain the Zeeman-splitting in m s"' and the magnetic field 
strength along the line of sight, which were determined from the 
cross-correlation between the RR and LL spectra. B|| was de- 
rived using a Zeeman-splitting coeflicient of 0.049 kms"' G"' 
(Vlemmings et al. i2006bj . 

4. Discussion 

4.1. CH^OH masers distribution 

The 6.7 GHz methanol maser emission likely arises when a 
methanol-rich cool gas (Tk = 30 - 100 K) at moderate densities 
(10-^ cm"-' < ne, < 10'' cm"-') is pumped by the infrared emis- 
sion from a layer of warm dust (Cragg et al. 120051 ). The warm 
dust is potentially associated with shocks that originate from the 
interaction of the VLA 1 jet and the surrounding medium. These 
are the shocks that give rise to the W75N(B) H2O masers that 
are found in elongated structures near VLA 1 (Torrelles et al. 
1997 ). The velocities of water masers are more extreme by about 
8 kms"' than those of the nearby methanol maser features Al- 
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Fig. 2. Total intensity and circular polarization spectrum for the 
brightest methanol maser feature of W75N(B). The thick red 
line in the bottom panel is best-fit total power derivative of a 
Gaussian fit of I spectrum to the circular polarization spectrum. 

A5, and the projected distance to VLA 1 of the H2O masers is 
shorter, suggesting that the methanol masers arise in the pre- 
shock region. 

As OH and methanol masers often occur in similar environ- 
ments, the methanol masers could be occurring in a structure 
related to the large OH maser ring around the three continuum 
sources (Hutawarakorn et al. |2002| l. The methanol masers dis- 
play a similar velocity gradient to the OH, with their velocity 
offset by only ~ 2 km s"' . However, no methanol emission is de- 
tected towards the southern part of the OH ring, making it more 
likely that the methanol maser structure is related to VLA 1 . 



4.2. The magnetic field 

Before we can relate the observed polarization to the magnetic 
field, we argue that the polarization vectors are perpendicular 
to the magnetic field and that the polarization is not aff'ected 
by non-magnetic phenomena. As in the case for H2O and SiO 
masers, the linear polarization of the methanol masers is par- 
allel or perpendicular to the magnetic field, depending on the 
angle between the magnetic field and the line of sight (ff). We 
solved the polarized radiative transfer using the equations from 
Nedoluha & Watson (1992) for the 6.7 GHz methanol maser 
transition. In Fig.[3]we show the relation between 6 and the frac- 
tional linear polarization. This relation depends on the emerging 
maser brightness temperature r^AO. We consider an upper limit 
of the maser beaming angle AQ < 10"' (Minier et al. ,2002j and 
the brightness temperature given from 



[K] 



S(y) 



[Jy] \[mas2] 



(1) 



where S(v) is the flux density, S the maser angular size, and 
^ a constant factor that includes all constant values, such as the 
Boltzmann constant, the wavelength, and the proportionality fac- 
tor obtained for a Gaussian shape by Burns et al. 1 19791 1. 



^= 13.63-10^ 



' Jy ' K. 



(2) 



Fig. 3. The angle 9 between the maser propagation direction and 
the magnetic field vs. the fractional linear polarization Pi for dif- 
ferent values of emerging maser brightness temperatures T^AQ. 
The dashed and dotted curves correspond to the brightness tem- 
perature limits for our observations (assuming AO. = 10"'). The 
horizontal solid curves show the fractional linear polarization of 
methanol maser features as reported in Table[T] The vertical solid 
line shows the most likely value of 6 for the W75N methanol 
maser observations. 



Since we find that all maser features are marginally resolved, we 
take a typical angular size value 1,-7 mas. Most of the maser 
features have T^AQ < lO' K sr, so we find that < Sent ~ 55° is 
ruled out across the entire methanol maser region. As shown 
in Goldreich et al. (I1973I I, this means that the linear polariza- 
tion vectors are perpendicular to the direction of the magnetic 
field. As indicated in Figl3] considering the observed brightness 
temperatures and assuming AQ = 10"', we estimate the aver- 
age angle between the magnetic field and the line of sight to be 
6 ~ 70°. For the non-paramagnetic molecules, however, there 
are a number of mechanisms that can create linear polarization 
and/or rotate the polarization angle. The incident infrared pho- 
tons on the masing region are anisotropically distributed, poten- 
tially producing high fractional linear polarization (P\), oriented 
tangentially to the incident radial photons. However, since the 
value of Zeeman frequency shift gQ. calculated for a magnetic 
field of 50 mG is greater than the rate of stimulated emission R 
(Vlemmings et al. I2009I I, the linear polarization vectors still di- 
rectly relate to the magnetic field. Consequently, in our case the 
magnetic field is oriented NE-SW with a position angle of about 
73° + 10°. The magnetic field is thus almost perfectly aligned 
with the large-scale molecular outflow, which is thought to be 
oriented close to the plane of the sky at a position angle of 66° 
(Hunter et al. '1994). Following our argument that the methanol 
masers are related to VLA 1 , the alignment of the magnetic field 
with the large-scale molecular outflows indicates that VLA 1 is 
likely the driving source of the outflow. 

The dust polarization of W75N has been observed at three 
wavelengths (450, 800, 1100 fim; Greaves & Holland fT998] l. 
These observations indicate a change in the magnetic field po- 
sition angle on the sky from 37° to -35° with increasing beam 
size. This change is attributed to the twisted magnetic field lines 
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around the massive protostar, a feature common to the models 
of low-mass star formation (e.g. Shu et al. 1995). In this case, 
the large beams are dominated by the dragged-in magnetic field 
lines, while the highest-resolution observations trace the compo- 
nents of the field that constrain the outflow (Greaves & Holland 
I1998I I. Our observations at mas-resolution are consistent with 
this model and suggest an outflow-driving mechanism similar to 
the one during low-mass star formation. 

The magnetic field derived from our EVN observations is be- 
tween 2 times larger than the one obtained using the Effelsberg 
telescope (Vlemmings 2008). This difference is expected, since 
the lower spatial resolution observations are often affected by 
spectral blending (e.g., Sarma et al. 2001) . Assuming 6 ~ 70°, 
the actual magnetic field strength in the methanol region is 
around 50 mG, which is similar to the region of the OH maser 
flare (Slysh & Migenes |2007] l. Crutcher (|1999) found an empir- 
ical relation between the magnetic field and the cloud density, 
B oc n^^, for number densities up to 10^ cm"-'. Later this relation 

was extended up to number densities of 10" cm"^ (Vlemmings 
12008 ). Hence, by applying this relation in the case of W75N, we 
find that the methanol masers occur in material ~50 times denser 
than OH, because Hutawarakorn et al. (120021) found a magnetic 
field ofVmG. 

The dynamical importance of the magnetic field can be quan- 
tified by considering the magnetic pressure and the dynamics 
pressure of the gas with velocity v and density p. Comparing the 
two pressures we can write 

ficrit = (3) 

which defines the critical magnetic field strength. If B < Bcrit, 
the dynamics of the region is regulated by the gas motion, and 
if B > B„h, the dynamics is dominated by the magnetic field. 
Considering a number density hh, ^ lO'^cm"'', above which 
the masers are occasionally quenched (Cragg et al. 120051) . and 
velocity of ~5 kms"' from internal maser motion, we obtain a 
Bait ^ 25 mG. As the magnetic field is larger than Bcrit over the 
whole 2000 AU area, it must play a dominant role in the dynam- 
ics of the W75N star-forming region. 
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5. Conclusions 

We observed the 6.7 GHz methanol masers in the massive 
star-forming region W75N with the EVN. We detected 10 
methanol maser sources, four of which are new detections. The 
linear polarization indicates that the magnetic field is aligned 
with the large-scale molecular bipolar outflow. This is consistent 
with the W75N(B) twisted magnetic field picture sketched by 
Greaves & Holland (1998). The detection of Zeeman-splitting 
indicates a magnetic field strength of the order of 50 mG. We 
suggest VLA 1 as the driving source of the large-scale molecular 
bipolar outflow of W75N(B). 
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